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Abstract 

The Magellanic Hi Stream (=s 2 x 10® Mq [c// 55 kpc]“) encircling the Galaxy at a distance d is arguably the 
most important tracer of what happens to gas accreting onto a disk galaxy. Recent observations reveal that the 
Stream’s mass is in fact dominated (3:1) by its ionised component. Here we revisit the origin of the mysterious 
Ha recombination emission observed along much of its length that is overly bright (~ 150 - 200 mR) for the 
known Galactic ultraviolet (UV) background («; 20 - 40 mR [iJ/55 kpc]^^). In an earlier model, we proposed 
that a slow shock cascade was operating along the Stream due to its interaction with the extended Galactic hot 
corona. We find that, for a smooth coronal density profile, this model can explain the bright Ha emission if the 
coronal density satisfies 2 x lO^'* < («/ cm“^) < 4 x 10“"^ at cl = 55 kpc. But in view of updated parameters for 
the Galactic halo and mounting evidence that most of the Stream must lie far beyond the Magellanic Clouds 
{d > 55 kpc), we revisit the shock cascade model in detail. At lower densities, the H i gas is broken down by 
the shock cascade but mostly mixes with the hot corona without significant recombination. At higher densities, 
the hot coronal mass (including the other baryonic components) exceeds the baryon budget of the Galaxy. If 
the Ha emission arises from the shock cascade, the upper limit on the smooth coronal density constrains the 
Stream’s mean distance to < 75 kpc. If, as some models indicate, the Stream is even further out, either the 
shock cascade is operating in a regime where the corona is substantially mass-loaded with recent gas debris, or 
an entirely different ionization mechanism is responsible. 


1. INTRODUCTION 

The Galaxy is surrounded by a vast amount of neutral gas 
in the form of high-velocity Hi clouds (HVC; Oort 1970). 
Eormally, these are neutral gas structures at a Galactic latitude 
\b\ > 30° having kinematic properties not consistent with the 
overall Galactic rotation (Wakker 2001). We now recognize 
that many of these make up the Magellanic Stream (MS; Di¬ 
eter 1971; Wannier & Wrixon 1972; Mathewson et al. 1974), 
roughly 2 x 10® Mq [c// 55 kpc]^ of gas (Eox et al. 2014) that 
has been stripped from the Magellanic Clouds (MCs), two 
dwarf galaxies in orbit around the Galaxy at a mean distance 
d^55 kpc (Walker 2012; Graczyk et al. 2014). 

The Magellanic Stream is ideal to study the environment of 
the Galaxy. Radio (Hi 21cm) surveys show that the Stream 
extends for 200° across the Southern Galactic Hemisphere 
(Nidever et al. 2010), and absorption line measurements to¬ 
wards distant quasars indicate a cross-section of roughly one 
quarter of the whole sky (Eox et al. 2014). Dynamical models 
(Besla et al. 2007; Guglielmo et al. 2014) constrained by ac¬ 
curate measurements of the proper motions of the MCs (Kalli- 
vayalil et al. 2013) agree that they move on a highly eccentric 
orbit, and that the MS spans a wide range in Galactocentric 
distance, from its source in the MCs system at roughly 55 kpc 
to 80 - 150 kpc above the South Galactic Pole (SGP) all the 
way to the tip of the tail. 

Given its relative proximity, the MS has been observed 
across the electromagnetic spectrum. Beyond H i, it has been 
detected in molecular (Richter et al. 2001) and ionized (Lu 
et al. 1994; Pox et al. 2005; Sembach et al. 2003) gas. A shad¬ 
owing experiment aimed at measuring the coronal soft X-ray 
emission discovered that the emission is enhanced in the direc¬ 
tion of the MS (Bregman et al. 2009). Recombination optical 
emission (Ha) was detected for the first time by Weiner & 
Williams (1996) and later confirmed by others (Reynolds et al. 
1998; Putman et al. 2004; Madsen 2012). But despite repeated 


attempts, to date, no stars have been discovered at any loca¬ 
tion along the Stream (e.g. Ostheimer et al. 1997). The data 
obtained by the recently completed Hubble Space Telescope 
(HST) Cosmic Origins Spectrograph (COS) UV absorption 
survey of the MS (Pox et al. 2013; Richter et al. 2013; Pox 
et al. 2014) indicates that the Stream is dominated by ionized 
gas, as was first proposed by Bland-Hawthorn et al. (2007). 
These data collectively suggest the existence of a strong in¬ 
teraction between the Stream gas and the hot halo (or corona; 
Spitzer 1956) of the Galaxy. 

Radiative hydrodynamic models (the ’shock cascade’; 
Bland-Hawthorn et al. 2007) indicate that the MS-halo inter¬ 
action may be strong enough to explain observed disruption of 
the Stream (Nidever et al. 2010) and its high ionization frac¬ 
tion. At the same time, the presence of coherent and strong 
enough magnetic fields (Putman et al. 1998) may stabilise 
the gas against severe ablation and provide thermal insulation 
to inhibit total evaporation of the neutral clouds (McClure- 
Griffiths et al. 2010). Whether such a shielding mechanism 
is operating all along the Stream is currently unknown. But 
without it, the Stream is likely to evaporate and mass-load the 
Galactic halo with a substantial amount of baryons (McCourt 
et al. 2015). 

Therefore, if we are to understand the complex environment 
of galaxies and how gas settles into galaxies (Heitsch & Put¬ 
man 2009), we need to explain first the observed properties of 
the Magellanic Stream within the framework of a multiphase 
hydrodynamical model. Elucidating the mechanism behind 
the bright spots of Ha emission observed along the Magel¬ 
lanic Stream has proven to be particularly challenging (see 
Bland-Hawthorn et al. 2007). To date, there have been two 
competing models (see also Konz et al. 2001): (i) the slow 
shock cascade discussed above; (ii) a new interpretation in¬ 
voking a powerful flare of UV radiation from the Galactic 
Centre (GC) powered by the accretion of material onto the 
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central black hole in Sgr A* (Bland-Hawthorn et al. 2013). 

The GC flare model is inspired by two important circum¬ 
stances: 1) The discovery of the y-ray emitting bubbles dis¬ 
covered by the Fermi satellite extending roughly 50° (10 kpc) 
from the GC (Su et al. 2010); 2) The observation that the 
brightest optical emission along the Stream is confined to 
a cone with half-angle 0i/2 ~ 25° roughly centred on the 
SGP (Madsen 2012). The GC flare model has found sup¬ 
port from the timescales and energy budget required to ion¬ 
ise the Stream, which are consistent with the results from jet- 
driven numerical models of the Fermi bubbles (Guo &. Math¬ 
ews 2012). More recently. Fox et al. (2014) have discovered 
that the ionisation levels over the SGP require an energetically 
harder ionising spectrum than elsewhere along the Stream, 
with the exception of a localised region near the LMC. 

The shock cascade model, on the other hand, explains the 
observation that the brightest Ha detections lie at the leading 
edges of the Hi clouds that make up the MS (Weiner et al. 
2002). However, this model may fail to produce the observed 
emission levels if the distance to the Stream at the SGP sig¬ 
nificantly exceeds the traditional view of if = 55 kpc, as in¬ 
dicated by most orbit calculations for the Magellanic Clouds 
over the past five years (Besla et al. 2012; Guglielmo et al. 
2014). The shock cascade model is strongly dependent on the 
density structure of the Galactic hot halo, and it assumes that 
the coronal density smoothly declines with Galactocentric dis¬ 
tance as cx r^^. But it now appears that both of these assump¬ 
tions may be false justifying our efforts to revisit the shock 
cascade model. 

The goal of this study is to investigate the strength of the re¬ 
combination (Ha) emission produced by the interaction of the 
MS gas with the Galactic corona exploring a range of Galacto¬ 
centric distances and different halo parameters, i.e., adopting 
different density profiles and temperatures of the gas sitting at 
rest in a fixed dark-matter (DM) potential. Note that through¬ 
out the paper we assume a flat, dark-energy- and matter (ba- 
ryonic and cold dark-matter; CDM) dominated Universe, and 
a cosmology defined by the set of parameters (relevant to this 
work) h - 0.7, G„, = 0.3, and Ga = 0.7. 

2. A MODEL OF THE GALACTIC HALO 


such a model given its solid physical foundation, with prop¬ 
erties that can be derived from first principles starting from a 
few basic assumptions (see e.g. Binney & Tremaine 2008). 

The scale-free potential, W, and the scale-free density, y, 
of an isotropic isothermal (DM) sphere are defined by (King 
1966): 


1 d I 2dW\ 
dx\ dx j 


9 yisoi^) 5 


y,so(x) = exp [W{x)] . 


( 1 ) 

( 2 ) 


where .x is a scale-free coordinate. 

This system of equations has no analytic solution, but it can 
be integrated numerically to values x <s: 1 to (nearly) arbitrary 
precision. Appropriate boundary conditions are, for example, 
the requirement that both the potential and the force vanish at 
the origin, i.e. VT(0) s 0 and W'{Q) = 0, respectively, which 
implies ^,so(0) = 1. Note that W{x) < 0 for x > 0. 

The connection between the scale-free quantities and their 
physical counterparts is given by 

r^r^x- p{r) = Pcy^soix) ; f{r) = (r^W{x) . (3) 


Here, pc, rc, and cr are the central density, the core radius, 
and the constant velocity dispersion, respectively. A solution 
corresponding to a particular physical system is obtained by 
flxing two of these three parameters (or any other two inde¬ 
pendent physical quantities of the system, for that matter); the 
third parameter is tied to the other two through the relation 
r^ - 9cp- 1 AnG Pc (King 1966). 

Alternatives to the isothermal sphere as a viable choice 
to describe a self-gravitating system of collisionless, DM 
particles, i.e. a DM halo include the Navarro, Frenk, & White 
(1996, NEW) model, 

i/NFw(0 = (^/^s)‘‘(l + '•/r9)^^ (4) 


where is a characteristic scale length; and the Einasto 
(1965) profile. 


ymAr) = exp ^ — 
a 



( 5 ) 


The density and temperature structure of the Galactic 
corona is largely determined by the underlying gravitational 
potential. The potential, in turn, is determined by the three 
main components of the Galaxy: the stellar bulge, the stel¬ 
lar and gaseous disc, and the DM halo. Given the mass of the 
bulge (~ 10*° Mo) and the disc (~ 10** Mg,; Kafle et al. 2014), 
and their size, these components are expected to dominate the 
Galactic potential only at r < 5 kpc and at r < 15 kpc, re¬ 
spectively. In other words, with exception of the inner ~ 5 kpc 
(Wegg et al. 2015), the DM halo dominates the Galactic poten¬ 
tial at all distances, with a similar contribution from the disc 
at scales comparable to its length. The distance to the nearest 
point of the MS is believed to be c/ ^ 55 kpc, which is the 
average distance to the LMC (« 50 kpc; Walker 2012), and 
the SMC (« 60 kpc; Graczyket al. 2014). Consequently, in 
modelling the interaction of the Magellanic Stream with the 
Galactic corona, it is safe to ignore the contributions to the 
potential from the bulge and the disc, and to focus instead on 
the DM halo only. 

2.1. The Galactic DM halo 

We model the DM halo of the Galaxy assuming it is well de¬ 
scribed by a single-component isothermal sphere. We opt for 


where a is a free parameter. Regardless of the model adopted, 
if the density profile p{r) is known, a general formalism can 
be applied which allows us to calculate scale-free quantities 
(q.v. Sternberg et al. 2002) on a case-by-case basis to describe 
a particular physical system. 

In order to specify our isothermal DM halo, and to compare 
its properties to the NEW and Einasto models, we proceed 
as follows. Based on the results of Kafle et al. (2014), who 
assume the dark halo of the Galaxy to be of NEW type, we 
fix the virial mass M„,> and the concentration x„,> s r„,> /r^. 
of the NEW dark halo to M„,> = 10*^ M©, and x„,> = 15, 
respectively. Our adopted value for the virial mass implies* 
r„ir ~ 260kpc and t;„,> s (GMcirlrdrY^^ ~ 130kms“*, re¬ 
spectively. Then we calculate the scale parameters for the 
isothermal and the Einasto DM halos by requiring that the 

* The viral mass and the virial radius are linked to one another through the 
relation 



We adopt a value for the cosmic mean matter density p,„ « 2.76 x 
g 337 Note that we define the virial radius at z = 0 (cf. 

Shull 2014). 
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Fig. 1.— Potential (left) and enclosed mass (right) of the Galactic DM halo. Each curve con'esponds to a different model, indicated by the legend in the left 
panel. The line of constant slope in the right panel shows the relation M„,>(r) given by equation (6). Note that the potential and the mass of all profiles have been 
matched at the virial radius of the Galaxy, r„,> k 260 kpc, indicated by the vertical line in both panels. See also Table 1. 


TABLE 1 

DM Halo Properties 



Isothermal 

Einasto^ 

NFW 

M„., [10*2 Mo] 

1.00 

1.00 

1 . 00 " 

tmr [kpc] 

259 

259 

259 

V„ir [km S“'] 

129 

129 

129 

^vir 

433 

13.8 

15 . 0 " 

t//„ir [10*^ cm2 g-2j 

- 1.11 

- 1.11 

-1.11 

r, [kpc]* 

0.60 

18.8 

17.3 

Vs [kms^‘]“ 

90.1 

82.4 

104 

MA(<r„,>)[10“> Mo]'* 

1.41 

3.03 

1.32 

Tom [105 K]“ 

5.81 

4.85 

7.73 


Note. — Fixed parameter values are shown in bold face; derived values, in 
normal font face. See text for details. “Gives the value of the ‘concentration’. 
The value x„,> = 15 is from Kafle et al. (2014). ^Corresponds to the scale 
radius for the Einasto and NFW profiles, and to the core radius for the iso- 
theimal sphere. “Con'esponds to the scale velocity for the Einasto and NFW 
profiles, and to the velocity dispersion cr for the isothermal sphere. '^Obtained 
from equation (10), assuming r = 1 and p = 0.59. “Obtained from equation 
(7), assuming p = 0.59. ^ We adopt u = 0.18 (Gao et al. 2008). 

virial mass and the physical potential at the virial radius in 
each case match the corresponding values for the NFW dark 
halo. Table 1 summarises our assumed values (in bold font 
face) and lists the derived values of the relevant scaling para¬ 
meters for each of the three halo models. The potential and 
corresponding mass for each model are shown in Figure 1. As 
can be seen, the relevant properties of an isothermal DM halo 
are very similar to the NFW and Einasto models. We note that 
the virial temperature is very similar across models, which is 
of relevance for the discussion in later sections. 

2.2. The Galactic corona 

It has long been known that the potential well of the Galaxy 
is filled with a diffuse, hot, gaseous component (the corona; 
Spitzer 1956). However, the origin of this gas, its thermody¬ 
namic state, its physical properties, its extension and hence 
its total mass, are still unknown. A recent attempt to con¬ 
strain the density structure of the corona implies a total gas 
mass of the Galaxy around ~ 10'° Mq within r„,> (Miller 
& Bregman 2013, 2015), which is consistent with other es¬ 


timates (e.g. Gatto et al. 2013). This important result relies 
on the assumptions that the corona is smooth, in collisional 
ionisation equilibrium, and isothermal, with a temperature of 
2 X 10® K. While the gaseous halo of the Galaxy is most 
likely not smooth nor strictly isothermal, the mean temper¬ 
ature of the gas inferred from its associated X-ray emission 
appear remarkably uniform across the sky (Henley & Shelton 
2013, 2014). 

Given this circumstance, we model the hot halo of the 
Galaxy as a single-phase, smooth, spherically symmetric com¬ 
ponent consisting of an ideal gas at a constant temperature T /„ 
in hydrostatic equilibrium with the DM potential i/r of an iso¬ 
thermal sphere. We further assume that the self-gravity of the 
gas is negligible, which is justified in the case of the Galaxy 
given that the inferred gas mass of the hot halo is on the order 
of 10“^ the mass of the DM halo (Suto et al. 1998). Under 
these assumptions, the total particle density of the corona is 
given by n(r) = noexp[i//(r)/a^]. Here, a^ = kThlpm,, is the 
isothermal sound speed, p is the mean molecular weight, 
is the atomic mass unit, and k is Boltzmann’s constant. 

The virial ‘temperature’ of the dark matter halo follows 
from the equivalence a^ = <t~, 


It worth emphasising that the above is is merely an equival¬ 
ent temperature, as p niu is not literally the DM particle mass 
(which is currently unknown). Table 1 lists the virial temper¬ 
atures of the different models. 

It is straightforward to show that under these assumptions 
the distribution of the hot gas in the potential is effectively 
governed by the thermal ratio (Cavaliere & Fusco-Femiano 
1976)2 

_ Tdm 2 / 2 /o\ 

T = — ^ o- la , (8) 

1 h 

such that 

n(r) = no exp [TW{r/rc)] , (9) 

^ Note that the designation of thermal ratio by the Greek letter j} is wide- 
spread in the literature. Here, we adopt its original designation. 
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Fig. 2. — Left: Density profile model (equation 9) adopting different DM halo models, all normalised to n{55 kpc) = 2 x 10“"^ cm“^. For a given halo, the curves 
correspond, from top to bottom, to values of t = 1.5, 1.0, and 0.5. Right: Gas mass enclosed within a given radius (equation 10) for each of the density profiles 
shown in the left panel. Note that the order of the curves is inverted with respect to the order on the left panel, i.e.. for a given halo model, the lowest curve 
corresponds to r = 1.5, etc. The vertical line indicates the virial radius of the Galaxy, « 260 kpc. Note that the enclosed mass for the Einasto model with r = 
0.5 is well above 10*^ M© and is not shown. 


The total gas mass within r for a particular value of r fol¬ 
lows from the integral of n over the appropriate volume. 


Mh(r) - Anfd 


X 


r 

n{x)x^dx. 


( 10 ) 


Clearly, the density field of gas in hydrostatic equilibrium 
with a fixed isothermal potential will be different for different 
gas temperatures. If the gaseous halo is ‘hotter’ than the DM 
particles, i.e. t < 1, then the density falls off more gently 
with radius. Values r > 1 on the other hand reflect the fact 
that the gas has cooled below the virial temperature of the DM 
halo, and is hence more concentrated. Large t values thus 
imply lower densities in the outer region of the halo, which 
in turn leads to weaker hydrodynamic interactions (at a fixed 
distance). 

Since the the virial temperature of our DM halo model is 
fixed, we control the temperature of the halo gas by vary¬ 
ing T. In the following we consider values of r in the 
range [0.5, 1.5] only, and choose the set of values r e 
(1.5, 1.0, 0.75, 0.5} as representative of this range. These val¬ 
ues imply gas temperatures T/, ~ 10® K (see Table 2), consist¬ 
ent with estimates of the temperature of the Galactic hot halo 
(Snowden et al. 2000). We ignore values of r < 0.5, since 
these yield overly shallow density profiles, which are incon¬ 
sistent with observations (see Figure 3). Similarly, we ignore 
values of T > 1.5 because these imply gas temperature which 
are too low compared to observations (see Section 3.2). 

We compare in Figure 2 the density profiles of gas at rest 
in different DM potentials, all scaled to a fiducial value n - 
2 X 10“^^ cm“^ at 55 kpc (see below and Section 3.2). For a 
given T, the gas density profile within a NFW DM potential 
(as defined by the respective parameters in Table 1) is steeper 
compared to the isothermal DM potential, but shallower than 
the Einasto DM halo. The difference in the density profile 
across the models reflects the difference in the potential in the 
radial range of interest. In contrast, the gas mass enclosed 
within r„,> for all three models and a given t are comparable 
(with exception of the mass for the Einasto model with t = 0.5 
which is well above 10*' Mq). Therefore, of all three models 
at a fixed t, the isothermal sphere halo leads to the highest 


gas density at any given distance beyond 55 kpc for roughly 
the same gas mass, and hence to the strongest hydrodynamic 
interaction in the outer halo. 

The different density profiles of gas sitting in an isothermal 
DM halo for our adopted values of t (equation 9), scaled to 
« = 2 X 10 cm“^ at 55 kpc, are shown in the top-left panel 
of Eigure 3. Eor comparison, we include there a set of val¬ 
ues of the halo density at various Galactocentric distances 
obtained from observations using a variety of methods (see 
Section 3.2). It is reassuring that all our models are fairly con¬ 
sistent with these measurements. 

The top-right panel of Eigure 3 displays the mass enclosed 
within a given radius (equation 10) for each of the mod¬ 
els shown in the top-left panel. As a consistency check, 
we compare the model masses to the upper limit on the 
Galactic baryion budget within r„,> set by the universal mean 
baryon-to-total mass ratio fb = Q.bl0.m- The most recent 
estimates of the baryon and cold dark matter mass densities 
= 0.02205 + 0.00028 and = 0.1199 + 0.0027 
(Planck Collaboration 2014) imply fb ~ 0.16. This value, to¬ 
gether with the total mass of the Galaxy {Mtot ~ 10'2 Mo) 
allow for a maximum gas mass of the Galaxy within r„,> be¬ 
low 10" Mq. Clearly, all our models result in masses within 
r„,> which are below this limit, indicated by the grey hatched 
area in the right panel of Eigure 3. Note that a density pro¬ 
file shallower than the r = 0.5 model, or a value of «(55 kpc) 
significantly higher than our fiducial value, would result in a 
gas mass largely inconsistent with these constraints. We in¬ 
clude in this figure the range of hot halo masses inferred from 
observations by Miller & Bregman (2015), comparable to the 
masses estimated by others (e.g Gatto et al. 2013), and the 
somewhat lower values inferred by Salem et al. (2015). Note 
that these estimates are all directly comparable to our model 
results since in all cases a smooth, monotonically decreasing 
density profile has been assumed. All our models, with excep¬ 
tion of the T = 0.5 model, predict masses within r„,> that are 
consistent with the inferred mass. 

Thus, the model r = 0.5 appears to be marginally consist¬ 
ent both with the mean density of the Galactic corona at large 
distances, and with the constraint on the total gas mass of 
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Fig. 3.— Left: Density profile model (equation 9) for different values of r. Also shown ai'e a set of halo density measurements at a range of Galactocentric 
distances. Right: Gas mass enclosed within a given radius (equation 10) for each of the density profiles shown in the left panel. The data points indicate the range 
of Galactic halo gas masses within r = 300 kpc infeiTed from observations. The grey hatched area on the top of the panel indicates the range of gas masses which 
are too large to be consistent with the (mean) cosmic baryon-to-total mass ratio ft ~ 0.16 (see text for details). Note that the top panels correspond to models that 
have been scaled to a fiducial value n(55 kpc) = 2 X 10^^ cm“^. The bottom panels show the corresponding results adopting fi(55 kpc) = 4 X 10^^ cm“^ instead. 
These are discussed in Section 4.4. 


TABLE 2 

Isothermal gaseous halo properties 


T 

Th [10« K] 

MA(A,r)[10“>Mo] 

0.5 

1.16 

3.36 

0.75 

0.77 

2.06 

1.0 

0.58 

1.41 

1.5 

0.39 

2.33 


the Galaxy. The model r = 1.5, although compatible with 
these constraints, appears too concentrated to be a plausible 
description of the Galaxy’s halo. In contrast, the models t = 
0.75 and r = 1 both display the best performance in terms of 
both the density profile and the gas mass enclosed within the 
virial radius of the Galaxy, although the former model yields 
a slightly hotter and more massive gas halo. In addition, these 
models reproduce by construction the relevant properties of 
the Galactic DM halo. Thus, we consider these models in 
particular provide a fully self-consistent and well founded de¬ 
scription of the Galactic hot halo, despite its idealized nature. 


3. NUMERICAL EXPERIMENT 

We simulate the passage of a stream of gas emulating the 
Magellanic Stream in its orbit through the Galactic hot halo 
by means of a 3-dimensional (3D) ‘wind-tunnel experiment’, 
expanding on the work by Bland-Hawthorn et al. (2007). In 
brief, we place a warm and essentially neutral, fractal gas 
cloud at a distance initially at rest with respect to the com¬ 
putational volume, and exposed it to a hot wind at a constant 
temperature Th (for a fixed t) and constant density n (for a 
fixed T and r^s) flowing with velocity fl/, under a fixed impact 
angle with respect to the gas cloud (see below). The warm 
gas is assumed to be initially in a state of pseudo-equilibrium 
with the hot gas, defined by the mean cloud-to-halo density 
ratio (or overdensity) rj = Pu,IPh^ and the mean cloud-to-halo 
pressure ratio f = PwjPh (see equation 1 1). 

3.1. Code 

We choose for our experiment the high-resolution, multi¬ 
phase, shock-capturing hydrodynamic grid-based code fyris 
ALPHA (Sutherland 2010), especially developed for astrophys- 
ical applications. The code solves the fluid dynamic equations 
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in one, two, and three dimensions as required. It has been 
shown to be fast, robust and accurate when compared to sim¬ 
ilar codes, and it performs well when subject to a standard 
suite of test cases as developed by Liska & Wendroff (2003). 
A unique feature of the fyris alpha code is that it includes non¬ 
equilibrium cooling through time-dependent ionisation calcu¬ 
lations. In addition, the code allows for the use of a variable 
equation of state (EoS) through a variable adiabatic index y 
and / or a variable mean molecular weight yU. These features 
are essential due to the large difference in the relevant time- 
scales which determine the physical state of multi-phase gas, 
as well as the large range of densities encountered in these 
type of simulations. 

3.2. Obsen’ational constraints, initial conditions, and 
set-up 

The proper motion of the Large Magellanic Cloud (EMC) 
has recently been measured using HST data, yielding an or¬ 
bital velocity Wmc = 321 +24 km s“* (Kallivayalil et al. 2013). 
We adopt the high-end value and set the speed of the hot wind 
to Vh - 350kms“*. Also, we adopt a value for the impact 
angle i? = 24°, such that the (shear) velocity of the hot wind 
is given by t)/, = (p;,cos??, p/,sin??, 0) = (320, 141,0)kms“*. 
Note that our results are fairly insensitive to the adopted value 
of §, as long as ?? » 0 and d <s: nl2, which is supported by the 
believe that the orbit of the Magellanic Stream is likely neither 
radial nor tangential with respect to the gaseous Galactic halo. 
We will assess the impact of this plausible, albeit arbitrary, 
choice on our results when dealing with virtual observations 
in Section 4 below. 

The Stream’s mean metallicity away from the MCs is now 
well constrained to Z » 0.1 Zq (Richter et al. 2013; Eox et al. 
2014). In contrast, the metallicity of the halo gas is still 
uncertain, although cosmological simulations (Rasmussen 
et al. 2009) and pulsar dispersion measures towards the EMC 
(Miller & Bregman 2015) both suggest that it is likely in the 
range Z ~ 0.1 Zq - 0.3 Zq far away from the disk. We choose 
a value for the metallicity of the halo of Z = 0.1 Zq, which 
is consistent with the mean value observed in external galax¬ 
ies similar to the MW (NGC 891; Hodges-Kluck & Bregman 
2013). 

Only gas clouds that are not overly dense and which have 
low pressure support with respect to the ambient medium 
will be disrupted in realistic timescales of ~ 100 Myr (Bland- 
Hawthorn 2009). We adopt ij - 100 and ^ = 0.1. With 
these parameters fixed, the initial temperature of the warm gas 
phase is set by the temperature of the hot gas phase through 

r. = (^)(f)A. (.1) 

where it should be noted that the mean molecular weight will 
be generally different in each phase. 

A key parameter of the models is the normalisation of the 
density profile at the canonical distance of the Stream above 
the SGP {d - 55 kpc). Although still uncertain, different 
lines of evidence indicate that it is likely in the range of 
10“^ cm“^ - 10“^ cm“^ at 20 kpc < r < 100 kpc. Eor ex¬ 
ample, Blitz &. Robishaw (2000) estimate a lower limit on the 
mean halo density of n « 2.4 x 10“^ cm“^ out to d < 250 kpc 
based on the assumption that the gas-poor dwarf spheroidals 
orbiting the Galaxy have been stripped from their gas by the 
ram-pressure exerted by the hot halo. Along the same line. 


TABLE 3 

Relevant simulation parameters / initial conditions 


Parameter 

Value 

Remarks 

(nj„ny,n-) 

(432,216,216) 

Grid dimensions 

{x,y,z) [kpc] 

(18,9,9) 

Physical dimensions 

6x [pc] 

42 

Spatial resolution (approximate) 

n [K] 

10*^ 

Halo gas temperature^ 

[K] 

10^ 

Initial Stream gas temperature^ 

M,(Hi)[Mo] 

10'' 

Initial Stream neutral gas mass 

V 

100 

Initial ratio of cloud : halo density 

( 

0.1 

Initial ratio of cloud : halo pressure 

n{55 kpc) [cm“^] 

2 X 10-* 

Total pai'ticle density at 55 kpc 

Zh [Zo] 

0.1 

Halo gas metallicity 

[Zq] 

0.1 

Stream’s metallicity 

X 

0.7154 

Hydrogen mass fraction 

Y 

0.2703 

Helium mass fraction 


24 

Impact angle 

Av [kms^*] 

200 

Velocity range of emission spectra 
with pixel size 6v = 2 km s~^. 

Note. — " Approximate. See Table 2. ^Approximate. The exact value 


will vary by small factors depending on Th (see equation 11). 


and combining observations with (2D) hydrodynamic simula¬ 
tions, Gatto et al. (2013) have inferred a range of halo dens¬ 
ities n M (1 - 4) X lO^^'cm'^ at 50 kpc < d < 100 kpc. 
Stanimirovic et al. (2002) have found that the gas clouds at 
the tail of the MS are likely in pressure equilibrium with the 
hot halo, and using this they have put an upper limit on the 
halo density of 10“^ cm“^ and 3 x 10“"^ cm“^ at a distance 
z = 15 kpc and z = 45 kpc from the Galactic plane, respect¬ 
ively. Anderson & Bregman (2010) infer a range for the mean 
halo density of n ^ (6 -10) x 10“^* out to the EMC {d ^ 5Q kpc) 
based on dispersion measures of EMC pulsars. However, the 
stripping of the LMC’s disc requires a somewhat lower value 
of ;?(48.2 + 3 kpc) = (1.1 ± 0.44) x 10“"^ cm“^ (Salem et al. 
2015). 

Based on these results, we adopt a fiducial value 
;?(55 kpc) = 2 X 10“^ cm“^, consistent with Bland-Hawthorn 
et al. (2007). As shown in Eigure 3, this choice leads to 
models for the Galactic corona that largely agree with the 
results from observations over a broad range in distances. 
In this respect, we consider both n{55 kpc) and r to be 
well constrained by observation. Note that the models are 
completely defined by the value of t, given that all the other 
parameters are either fixed or they depend on t (Table 3). 

We run all simulations in a rectangular box of comov¬ 
ing size 18 X 9 X 9kpc^, using a fixed grid composed of 
432 X 216 X 216 cells. These settings imply a spatial resol¬ 
ution of 6x - (9/216) kpc ai 42 pc. The fragment of gas 
representing the Magellanic Stream is initially constrained to 
a cylinder 18 kpc in length and 2 kpc in diameter, and whose 
axis of symmetry runs parallel to the x axis of the coordin¬ 
ate system defined by the box. The impact angle is defined 
with respect to the x-axis of this cylinder. In this setup, the 
x-axis coincides with the Magellanic longitude, Im, whereas 
any of ?/ or z run along the Magellanic latitude, (Wak- 
ker 2001; Nidever et al. 2008). The simulated Magellanic 
Stream consists of an H i gas distribution initially at temper¬ 
ature Tu, ~ 10^ K (see equation 11); a mean initial hydro¬ 
gen particle density n ~ 10“^ cm“^; and a total neutral gas 
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mass ~ 10^ Mq. The initial warm gas density field 

corresponds to the density of a fractal medium described by 
a Kolmogorov turbulent power spectrum P{k) oc with 

a minimum wavenumber k,„i„ — 8 (relative to the grid) cor¬ 
responding to a spatial scale of 2.25 kpc (q.v. Sutherland & 
Bicknell 2007), comparable to the typical size of clouds in 
the Stream. The H i gas cloud is assumed to be at a fixed dis¬ 
tance Tms from the Galactic Centre in the direction of the SGP. 
The temperature Th and the density n of the hot wind are set 
according to the value of t and the distance r^s as given by 
equations (7) and (9), respectively. 

For each r 6 {0.5, 0.75, 1.0, 1.5), we consider a set of 
Galactocentric distances e {55, 75, 100, 125, 150} kpc 
which together span the range of plausible orbits of the MS 
above the SGP (Guglielmo et al. 2014). This yields 20 
models. Each model is run for a total (simulation) time of 
320 Myr, starting from fs„„ = 0 Myr, assuming free bound¬ 
ary conditions. The simulation output for a given set of val¬ 
ues {t, r„s, tsim] - in steps of Atsm - 10 Myr - consists of a 
series of datacubes containing information about the H i- and 
Hii densities, n hi and riuu, respectively; the gas temperature 
T, and the gas velocity v = {Vx, Vy, v^). Using this informa¬ 
tion, we compute the Ha emission and the H i column density 
of each cell, and from these the Ha surface brightness and 
total HI column density along the sightline for each snapshot. 


3.3. Emission line spectra 

We compare the Ha emission of the gas in our simulations 
to Fabry-Perot Ha observations along the Magellanic Stream 
(Bland-Hawthorn et al. 2013, and references therein), and 
complement these with results on the associated H i column 
density measurements. 

Given the one-to-one correspondence between the H 1 21 cm 
emission (i.e., brightness temperature) Tb and the Hi column 
density of a parcel of optically thin gas (e.g. Dickey & Lock- 
man 1990), we use the total H i column density along the sight¬ 
line (Nhi) as a proxy for the corresponding Hi 21cm emis¬ 
sion, i.e. we define the H i 21cm intensity to be /hi = Nhi. 

The Ha emission is computed using 


PHa 


PHa 


{shock) 


+ Pna 


{phot) 


( 12 ) 


The first term accounts for the ionisation that results from 
slow shocks produced by the collision of the trailing cloud 
gas with the leading gas ablated by the interaction with the 
hot halo, and is given by 

where is the effective Ha recombination coefficient 

(equation Al), Kr ^ 1.67 x 10 cm^ s mR,^ and the factor 
[1 + {YjAX)} accounts for the conversion of electron density 
to HII particle density. Adopting a hydrogen and helium mass 
fractions X - 0.7154 and Y - 0.2703, respectively (Asplund 
et al. 2009, solar bulk composition),"^ assuming the gas is fully 
ionised («„ » nu,,), it follows that [1 + {YjAX)] » 1.09.^ 


The second term in equation (12) accounts for the ionising 
effect of the cosmic ultraviolet background radiation (UVB). 
The Ha emission along the sightline of gas in photoionisation 
equilibrium with the UVB radiation field is 

[faa hhi) ds. (14) 

where fna{T) gives the fraction of recombinations that pro¬ 
duce an Ha photon, and /Hff(10"^K) 0.45 (equation B2). 

The HI recombination rate, Thi, is related to the total ionising 
photon flux O, through 

O,- = 1.59xl0Vhotoncm-2s-‘ ■ (15) 

We adopt Fhi = 10 s“' (appropriate for z = 0; Weymann 

et al. 2001) and y = 1.8 (Shull et al. 1999), corresponding 
to an ionising flux <5, ~ 10"^ photons cm“^ s *. If we used 
instead the most recent estimate Thi = 4.6 x 10“^"^ s“' (Shull 
et al. 2015), the flux would be lower by roughly a factor 2. 

To mimic radiation transfer effects, we limit the depth 
(along the sightline in any direction) of the gas ionised by 
the UVB to a maximum value defined by the condition that 
the column recombination («e «hii Q'b) equals the incident ion¬ 
ising photon flux (O,). This condition is equivalent to restrict¬ 
ing the ionising effect of the UVB to a column of neutral gas 
~ 10*^ cm“^ (see Appendix B). The effect of the cosmic UVB 
is to produce an ionisation skin around the cloud featuring an 
Ha surface brightness at a level of roughly 5 mR. Again, if we 
used instead Thi = 4.6 x 10 s ', this value would decrease 
to roughly 2 mR. It is important to mention that this approach 
is not entirely self-consistent with our simulations because the 
ionising effect of the UVB is not included at runtime, and be¬ 
cause it assumes photoionisation equilibrium. Also, we ig¬ 
nore for the moment the contribution of the Galactic ionising 
field, which would produce an additional mean Ha signal of 
21 ^ {dl55 kpc)“^ mR (^ 2; Bland-Hawthorn et al. 2013). 

To allow for a faithful comparison with observations, we 
map the simulation data onto observed space by projecting 
the simulation volume along a given axis, so as to mimic the 
projection of the observed emission along the Stream onto the 
plane of the sky. We choose, for convenience, an axis parallel 
to one side of the simulation box, and perpendicular to the 
Stream’s main axis, i.e. the i/-axis. However, we will address 
the potential bias introduced by this choice by comparing the 
results obtained by projecting along all three orthogonal axes, 
X, y and z- 

Each sightline across the projected datacube thus corres¬ 
ponds to a pencil-beam spectrum. Ha and Hi 21cm pencil- 
beam spectra along the chosen projection axis are construc¬ 
ted by computing the intensity for each cell in the simulation 
volume. The Ha intensity and the Hi 21cm emission intens¬ 
ity observed at velocity y,,, (in the rest-frame of the Hi gas 
at tsim = 0) of a parcel of gas at cell n s (/, j, k) with bulk 
velocity Va are, respectively, 

~ tt) — ^Ho'(tt) ^n) 5 

7Hi(^m “ n) = N Hi(tt) <pn{^m ~ ^n) 5 (16) 


^ 1 milli-Rayleigh (mR) corresponds to 10^/4;r photons s“* sr“* 

(Baker & Romick 1976), or 2.41 X 10“'* erg s“* sr“' at Ho'. 

For our adopted metallicity of 10 percent the solar value, Zq = 0.0142, 
the contribution of heavy elements to the electron density can he neglected. 

^ We assume that helium is only singly ionised, given that the ionisation 
energy of Hell is comparatively high (E x 54.4 eV; Kramida et al. 2014). 


where the normalised line profile (piv) - exp with 

bj = (2kTJm,fil\ 

We adopt a spectral range in terms of velocity of Av - 
200 km s ' which corresponds to the spectral range provided 
by the WHAM spectrometer. The velocity scale is given with 
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Fig. 4. — Virtual observation of the Ha (left) and Hi 21cm (right) emission along the Stream in model t = 1 and = 55 kpc at f,,-,,, = 170 Myr, projected 
along the ij axis. Top: Emission at the intrinsic (high) resolution of the simulation. The arrow on the top-left comer indicates the projected Stream’s velocity. 
Middle: The result of spatially smoothing the emission map using a circular Gaussian kernel of size 0.5° / 1° (FHWM / full width), as indicated by the circle on 
the bottom-right corner of the top panel. Bottom: Emission at different pointings obtained using a beam with an angular diameter of 1° on the sky. The position 
of each individual pointing is dehned by each cell of the grid overplotted on the middle panel. The grayscale bar to the right of each set of panels indicates the 
Ha intensity in milli-Rayleigh (left) and the H i column density in cm“^ (right). The full evolution of the Ha surface brightness and the H i column density are 
available at http://www.physics .usyd.edu.au/-tepper/proj_ms .html. 


respect of the initial rest-frame of the H i gas, such that emis¬ 
sion spectrum spans the range [-100 ,-hIOO] kms“*, with a 
pixel size 5v - 2kms“*. 

At each beam position (or ‘pointing’), we compute the 
intensity-weighted average of all cells within the beam, res¬ 
ulting in a single spectrum per pointing. We adopt a low 
resolution beam of 1° diameter on the sky, identical to the 
resolution provided by the WHAM spectrometer. 

We choose a rectangular (rather than a circular) beam, 
which allows for a full coverage of the projected image, and 
which greatly simplifies the scanning procedure.® Given 
that a beam with a diameter of 20 subtends a solid angle 
^beam ~ 7r(0/2)^ (provided that 0 <sc 1), the linear size of 
a square subtending a solid angle Q.beam at a distance r„s relat¬ 
ive to the angular dimension of a single cell in our simulation, 
dfiy (djc/ r„s)^, is roughly 61 - 0. The 

emission within a beam pointing at each velocity v,„, i.e. the 
beam spectrum, is then 

[<5/1 [<5/] 

* A alternative arrangement consisting of a tightly packed array of circular 
windows (which conveys equal weight to every pixel within the beam; see 
e.g., Haffner et al. 2003, their Figure 3) would yield essentially the same 
results. 


where the sum extends over all cells within iibeam, and X e 
(Hi, Hii). Here, the notation [k] indicates the largest odd in¬ 
teger smaller than or equal to k. 

The average of the emission within each beam pointing can 
effectively be obtained by overlaying a rectangular grid on the 
projected image with a cell size equal the solid angle subten¬ 
ded by a circular beam of diameter 20 at that distance, and 
computing the arithmetic mean within each new cell. We 
choose the origin of the matrix to be shifted by half a beam 
size in each direction (xy) to avoid the uncertainties associated 
with the simulation volume’s boundaries. 

It is important to emphasise that the pixelation of the sim¬ 
ulation volume and the position of the beam pointings with 
respect to the projected datacube are rather arbitrary. Also, 
the latter is also generally different for each adopted distance 
r„s. This results from the fact that for a computational volume 
with fixed comoving size and fixed grid, and a beam of fixed 
angular size, an increasing fraction of the gas that represents 
the Stream will be sampled by the beam with increasing dis¬ 
tance. Because of this, and also to avoid a bias in the result¬ 
ing emission introduced by potentially bright features induced 
by chance alignments (rather than due to intrinsically bright 
gas blobs), prior to computing the spatial average within each 
beam pointing we smooth each 2D spatial slice at each ve¬ 
locity bin using a circular Gaussian kernel with a full width 
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at half maximum of half the beam size and a total width of 
(i.e., truncated at) the size of the beam. According to the ap¬ 
proach described above, the projected intrinsic Ha emission 
map from our simulation is smoothed at each given distance 
using a Gaussian kernel with FWHM of [5//2] pixel and a 
total width of {Sl\ pixel. 

Finally, in order to take into account the instrumental line 
broadening, we convolve each spatially averaged spectrum 
using a Gaussian kernel^ with a FWHM = 12kms ’, which 
roughly corresponds to the resolution of the WHAM spec¬ 
trograph (Reynolds et al. 1998). The Ha surface brightness 
and Hi 21cm emission maps are simply obtained from the 3D 
spectral datacube by simply integrating each spectrum along 
the velocity coordinate (equations 16). 

A selected example illustrating the result of the above pro¬ 
cedure is shown in Figure 4 (see also Figure 6). As we show 
below (see Section 4.3; Figure 7), the Hii to Hi mass ratio 
in this snapshot roughly matches to the corresponding ratio 
observed in the Magellanic Stream (~ 3; Fox et al. 2014). 
The effect of the ionising cosmic UVB is apparent: all the 
Stream gas is lit up and emitting at level of ~ 5 mR. But there 
are brighter spots which are a consequence the shock cascade, 
whereby the trailing clouds collide with the material ablated 
by hydrodynamic instabilities from the leading gas, thus be¬ 
ing shock ionised (q.v. Bland-Hawthorn et al. 2007). Interest¬ 
ingly, while the bright Ha spots seem to closely track the high 
N Hi parcels of gas, the converse is not true, with high N hi ap¬ 
pearing with no correspondent strong Ha emission. However, 
these differences become less apparent, although they remain, 
as a result of the beam smearing. The most dramatic effect of 
the latter is the dilution of the Ha and H i signals with respect 
to the brightest levels seen at the (intrinsic) resolution of the 
simulation by a significant factor. Only the brightest spots in 
Ha would observable with a Fabry-Perot interferometer for 
reasonable integration times, leading to a significant fraction 
of the ionised gas mass falling below the detection threshold. 

4. RESULTS 
4.1. Gas emission 

We follow the evolution of the gas emission over a period of 
320 Myr in all 20 models. To illustrate these results, we adopt 
the model t = 1 as our standard model. The top panels of 
Figure 5 shows the evolution of the Ha / H i emission at four 
representative distances r„s = 55 kpc, 75 kpc, 100 kpc, and 
150 kpc. The middle panels show the corresponding result for 
all models r = 0.5, 0.75, 1.0, and r = 1.5, at a fixed distance 
r„s = 55 kpc. Note that the comparison of the results for a 
fixed T at different distances allows to assess the effect of the 
density on the Ha / H i emission for a fixed temperature, while 
the comparison of models with different t at a fixed distance 
of 55 kpc (and thus a fixed density) helps us explore the effect 
of the temperature. 

We find that the gas ionises quickly (~ 50 Myr), and after ~ 
300 Myr, the H i column density has decreased uniformly both 
with time, dropping by nearly an order of magnitude. This can 
be understood as a consequence of the increasing ionisation 
of the gas due to the interaction with the hot halo gas. The 
detected column density of the gas farther out is also lower 
with respect to the intrinsic gas density, which is an effect of 

^ Note that the LSF of the WHAM spectrometer is only poorly approxim¬ 
ated by a Gaussian (Tufte 1997). However, given the typical Ho- / Hi 21cm 
line widths (20 km s“* - 40 km s“*), this approximation hai'dly affects our 
results. 


the beam dilution. Note that a beam of 1° diameter samples a 
region of roughly twice (three times) the size at = 100 kpc 
(150 kpc) with respect to r„s = 55 kpc. 

Similarly, the highest Ha emission - this is, the maximum 
value at each given time - comes from the gas which is closest. 
In contrast to the behaviour of the H i density, in this case the 
effect is governed by the increase of halo gas density with 
distance, and the corresponding strength of the hydrodynamic 
interaction leading to the shock cascade. What is surprising is 
that even at the lowest Galactocentric distance of r„s = 55 kpc, 
the emission never exceeds ~ 40 mR; and it barely reaches 10 
mR at 75 kpc. At even larger distances, r„s ^ 100 kpc, the 
emission is dominated by the recombination of the gas ion¬ 
ized by the cosmic UV background. Note that often we do 
find in our simulation pixels with //ho > 30 mR, and occa¬ 
sionally on the order of ~ 100 mR, but their strong signal is 
washed out as a result of the beam smearing (Figure 4; see 
also Section 4.4). The similarity in Ha emission across mod¬ 
els with different r shows that these results are insensitive to 
variations in the halo gas temperature by factors of a few. This 
indicates that the Ha emission is dominated by the gas ionised 
through cloud-cloud collisions that trigger the shock cascade, 
rather than the gas ionised by the interaction with the hot halo. 
The insignificance of the halo gas temperature together with 
the fact that the virial temperature across the DM halo mod¬ 
els presented previously is very similar (Table 1) makes the 
choice of the DM halo model irrelevant, as long as the corres¬ 
ponding gas density profiles are comparable. 

We see that the maximum level of Ha emission at any 
reasonable distance is comparable to, or even less than, the 
emission induced by the Galactic ionising starlight ~ 20 - 
40 mR [c//55 kpc] “. This is the reason for us to ignore this 
component in our models, since it would otherwise outshine 
the emission produced by the shock cascade. Taking the con¬ 
tribution of the Galactic UV into account would elevate the 
emission at 55 kpc to ~ 70 mR, and to ~ 30 mR at 75 kpc, 
which are significantly lower than the Stream’s emission ob¬ 
served over the SGP and at the tail of the Stream at Im ~ 260°, 
respectively (Bland-Hawthorn et al. 2013). It is however un¬ 
clear at this point how the ionising effect of the Galactic UV 
included self-consistently at runtime would affect these limits. 

Given that observations performed with the WHAM instru¬ 
ment typically reach a sensitivity of > 30 mR, the Stream gas 
in our models would be essentially undetectable (ignoring for 
the moment the contribution which results from the Galactic 
ionising field). In contrast, and considering that the sensitivity 
of e.g the GASS survey is roughly Nhi = 1.6x 10** cm“^, the 
gas at distances if < 100 kpc would be bright in H 1 21 cm, and 
marginally detectable at ~ 150 kpc, even after 300 Myr. 

We find that much of the gas dislodged from the main 
body of our model Stream is low density material that mixes 
rapidly with the halo gas, thereby being heated (and thus 
ionised) to temperatures well above 10^ K, which are on 
the order of the temperature expected for turbulent mixing 
(Begelman & Fabian 1990). At these temperatures, the Ha 
emissivity drops by nearly two orders of magnitude with 
respect to its value 10"* K (see equation Al), and the ionised 
gas becomes thus practically invisible in Ha. Only the 
gas ionised by cloud-cloud collisions remains at relatively 
low temperatures (~ 10"* K), and recombines quickly, thus 
providing the strongest Ha signal. However, the fraction 
of warm ionised gas is very low overall, and thus is the 
corresponding Ha signal. Hence, only a mechanism such as 
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Fig. 5. — Evolution of the Hi column density (left) and the Her surface brightness (right). The top panels correspond to the results for model r = 1, and a 
range in Galactocentric distance r^s (indicated by the legend). The middle panels show the result at a fixed Galactocentric distance = 55 kpc, for different r 
(indicated by the legend). Each data point in the left (right) panels corresponds to the 90 (100) percentile of the distribution of intensities at that particular time 
and distance. The four top panels correspond to the results obtained from projecting the simulation cube along the y-axis. The two bottom panels illustrate the 
difference between projections along all three orthogonal axes for model r = 1 at 55 kpc. 
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slow shocks which is able to ionise a significant fraction of 
the gas without increasing its temperature far above 10"^ K 
will lead to significant Ha emission. 

As the reader may recall, all the above results correspond to 
virtual observations where the simulation cube has been pro¬ 
jected along the j/-axis. Given that this choice is somewhat ar¬ 
bitrary, we calculate the corresponding results for projections 
along all three orthogonal axes for model r = 1 at 55 kpc. The 
outcome of this exercise is summarized in the bottom panels 
of Figure 5. We have checked that the results are essentially 
the same for all other models. The Hi column density is con¬ 
sistently highest when observed along the x-axis, given that 
this axis coincides with the axis of symmetry of the initially 
cylindrical gas configuration in our setup, and the sightline 
traverses a larger path across the Hi cloud. The projections 
along the y and z axes yield nearly identical results, as expec¬ 
ted from the symmetry of the initial cloud structure. 

In our experiment, the projection along the x-axis is equival¬ 
ent to observe the Stream ’face-on’. Hence, one would naively 
expect that the Ha emission should be highest when project¬ 
ing along this direction. But surprisingly, the Ha intensity 
in our models is very similar regardless of the projection, be¬ 
ing only slightly stronger when viewing the gas cloud face- 
on. This implies that the choice of projection axis to compute 
virtual observations is essentially irrelevant. Moreover, since 
different projections effectively imply significant variations in 
the impact angle, our choice of a particular value for i? turns 
to be irrelevant as well, as far as the Ha intensities are con¬ 
cerned. 

4.2. Gas kinematics 

The kinematics of the warm neutral and ionised gas phases 
as traced by Ha and H i emission provide a insight into the 
mechanism ionising the Stream. We explore this using the 
spectra of our model identified by r = 1 and r„s = 55 kpc at 
tsim - 170 Myr. Note that the results, with exception of the 
Ha emission strength, are virtually identical for all other mod¬ 
els. An example of a strong Ha and Hi 21cm emission lines 
is shown in the top panel of Figure 6. These correspond to the 
pointing with the brightest Ha emission shown in the bottom 
panel of Figure 4. The integrated strength of the line is indic¬ 
ated in each case in the top-left corner. For reference, the typ¬ 
ical sensitivity of WHAM (30 mR) translates into a spectral 
sensitivity of 0.5 mR/ kms“\ assuming a typical line width 
of 30 km s“'. Similarly, the sensitivity of the GASS survey 
( Nhi = 1.6 X 10** cm“^) corresponds to a spectral sensitiv¬ 
ity of 10*’ cm“^/ km s”*. Thus, while the Ha emission in this 
case is marginally above the WHAM detection threshold, the 
corresponding Hi 21cm signal would be comfortably detec¬ 
ted in a survey similar to GASS. We find typical line widths 
of 20 - 30kms * (FWHM), which are consistent with obser¬ 
vations (Putman et al. 2003). 

We study the difference in the kinematics of the warm 
neutral and ionised gas by comparing the sightline velocity 
centroid of the Ha emission to the Hi 21cm emission, dis¬ 
tinguishing between pointings ‘on’ and ‘off’ the Hi clouds. 
In this context, ‘on’ (‘off’) means that the Hi 21cm emis¬ 
sion is above (below) the GASS detection limit Nhi = 
1.6 X 10** cm“^. Note that the line centroid corresponds to 
the intensity-weighted mean velocity. In addition, we flag 
those pointings where the Ha emission is above the level ex¬ 
pected from ionisation by the UBV (Figure 6, bottom panel). 


We do not find a significant difference between the sightline 
velocities of the warm neutral and ionised gas phases; their 
respective velocity centroids agree within +5kms“*, as ob¬ 
served (Putman et al. 2003; Barger et al. 2015). There is a 
tendency for the Ha lines to have higher slightly higher ve¬ 
locity centroids. This arises from the higher line asymmetry 
resulting from a more extended emission along the sightline. 
It is interesting that we barely find any Ha emission above 5 
mR detached from Hi 21cm emission. Also, the velocity of 
this strong Ha emission is generally low (y < 20 km s '), with 
a tendency for the strongest emission to have the lowest velo¬ 
cities (not shown), indicating that the strong Ha emission is 
physically associated to the Hi gas. This coincidence in both 
velocity and physical space of the Ha/ H i signal is - recall the 
weak dependence of the Ha emission on the halo temperature 
- another characteristic signature of the shock cascade. 

4.3. Gas ionisation timescales 

It has recently been inferred that the mass of the ionised gas 
kinematically associated to the Magellanic Stream is roughly 
3 times larger than its Hi mass (Fox et al. 2014). Here, we 
briefly explore the evolution of the ionised and warm-neutral 
gas mass fractions, using the H ii and H i masses in as a proxy 
for the ionised and warm-neutral gas phases, respectively. 
Again, we focus on the results obtained from our models char¬ 
acterised by T = 0.50, 0.75, 1.00 and r„s = 55 kpc. Note 
that the results are qualitatively the same for all other models. 
Figure 7 shows the evolution of the individual mass fractions 
fx = MxIMfoi, where X e { H i, H ii} and = Mh, + Mh„, 
as well as the evolution of the H ii to H i mass ratio. We in¬ 
clude for reference the inferred value of the ratio Hii ; Hi 
= 3. We find that the gas evolves on a typical timescale of 
100 Myr, which is consistent with the time estimate which 
results from assuming that the cloud disrupt by the action of 
Kelvin-Helmholtz instabilities and our adopted (initial) value 
of Tj. After this period, the ionisation effect of the halo-cloud 
and cloud-cloud interaction leads to a reduction of the warm 
neutral gas mass by half. The Hii ; Hi mass ratio increases 
rapidly with time, and the inferred Hii; Hi mass ratio of 3 is 
reached after ~ 170 Myr. 

Since our simulations assume free boundary conditions, it 
is difficult to quantify with precision how much of the neutral 
gas is ionised and how much simply escapes the simulation 
volume. Nonetheless, we estimate that only a negligible 
fraction of the warm-neutral gas is lost by fj„„ < 270 Myr. 
Therefore, the ionisation of the gas due to interactions with 
the hot halo gas, and due to cloud-cloud interactions lead 
to a strong evolution of the mass fractions in the neutral 
and ionised phases. The relatively short survival timescale 
implies the requirement for a continuous replenishment of 
gas from the MC to the Stream (Bland-Hawthorn et al. 2007). 

Note that, for the ease of discussion, in the following we 
shall refer to the set of models discussed in the last sections, 
defined by n(55 kpc) = 2 x 10“^* cm“*, as the ‘standard’ model 
set. 

4.4. Conservative departures from the standard models 

Given the results of the last sections, our standard set of 
models appear to indicate that the shock cascade fails to 
produce the mean level of Ha emission (~ 160 mR; Bland- 
Hawthorn et al. 2013) observed along the Magellanic Stream. 
But there are two factors that deserve closer consideration. 
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Velocity (km/s) 



Fig. 6. — Gas kinematics in model r = 1 and r^s = 55 kpc at = 170 Myr. Left: Emission spectra observed at the pointing with the brightest Htr emission 
shown in the bottom left panel of Figure 4. We show for comparison the intrinsic line profile and the line profile convolved with a Gaussian LSF with a 
FWFIM = 12 kms^*. Right: Line-of-sight velocity centroid of the Flo- intensity detected 'on' and ’off’ FIi clouds (see text for details). The solid diagonal lines 
indicate the identity line and a ±5 kms“* range around this line. 



Time (Myr) 

Fig. 7.— Evolution of the FIi and Fin mass fractions in the r = 1 model 
Stream at r^s = 55 kpc, projected along the y-axis. The top panel shows the 
evolution of the individual mass fractions Mx/^wt, where X € { FIi, FIii| 
and = Mhi| + Mnn- The bottom panel shows the evolution of the Fin 
to FI I mass ratio. The dashed horizontal line indicates the ratio H ii : FI i = 3 
infen'ed from observations (Fox et al. 2014). The results are essentially the 
same for all other models, and are therefore not shown. 

First, a halo density at any given distance within r„,> higher 
than implied by our standard models could enhance the on¬ 
set and development of hydrodynamic instabilities (Kelvin- 
Helmholtz), thus promoting the shock cascade and the result¬ 
ing Her emission. Secondly, an increase in the beam resolution 
would certainly diminish the smearing effect on bright spots 
which have characteristic sizes significantly smaller than the 
beam. Consider, for instance, that a beam with a diameter of 
10’ samples a region which is nearly ten times smaller than 
a 1° beam at r^s = 55 kpc. Indeed, the brightest Ha observa¬ 
tions along the Stream have been obtained with spectroscopy 
over smaller apertures (3’-10’) than the WHAM survey (e.g. 
Putman et al. 2003). 

Therefore, it is important to extend the parameter space 
of our study, in terms of both the halo density and the ad¬ 
opted beam size. We now increase the normalisation of the 


halo gas density at 55 kpc by a factor of 2, i.e. n(55 kpc) = 
4 X 10“^ cm“^, but with initial conditions and set up which 
are otherwise identical in every aspect to the standard mod¬ 
els. These shall be referred to as the ’extended’ models. In 
addition to increasing the density, we produce a new set of vir¬ 
tual observations for the standard models, adopting a smaller 
beam size (i.e. a higher resolution) with a diameter of 0 = 10' 
(rather than 1 °), and two sets of virtual observations for the ex¬ 
tended models, adopting either a low (1°) or a high (10’) res¬ 
olution beam. We shall refer to these as the ‘low-resolution’ 
and ‘high-resolution’ models, respectively, keeping in mind 
that is not the actual hydrodynamical model, but the virtual 
observation, to which the resolution refers. 

Note that the density profiles implied by the extended mod¬ 
els are fairly consistent with observations (Figure 3, bottom- 
left panel), although there is no model which agrees with the 
data over the full range in distance. However, the enclosed 
gas mass in the halo within r„,> resulting from each of these 
density profiles is larger than the the mean range of masses 
inferred from observations, and - with exception of the t — I 
extended model - they are all inconsistent with the mass limit 
imposed by the universal baryon-to-total-mass fraction (Fig¬ 
ure 3, bottom-right panel). Hence, all the extended models 
but the T = 1 model, may be deemed ’unphysicaF. Nonethe¬ 
less, it is still of interest to explore the Ha intensity in these 
type of models, as will be discussed later. 

Each of the models in either the standard or the extended 
set is run for a total simulation time of tsim - 320 Myr, and 
a virtual observation of the Ha intensity of the gas at the ap¬ 
propriate resolution is produced every Af,™ = 10 Myr. There¬ 
fore, for each model and at each time step we obtain a whole 
new distribution of Ha intensities. In order to deal with the 
overwhelming amount of information, and to make a mean¬ 
ingful comparison between models, we opt for the following 
approach: Since the total time lapse fj™, and the choice of 
output time step are somewhat arbitrary, for each model we 
single out the snapshot at which the maximum Ha intensity 
anywhere in the gas (i.e. at any beam pointing) is largest. This 
is further justified by the fact that the maximum Ha emission 
does not evolve strongly with time (Figure 5). Note that the 
snapshot thus selected will in general be different for each 
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Fig. 8.— Overview of model results. In any panel, a model is identified by the set of parameters {«(55 kpc), B, r, rMs|. The corresponding data point consists of 
a symbol (• or x), and two numerical values. The value to the top-right of a given symbol (in parentheses) indicates the maximum Hn emission (mR) observed in 
the full time lapse of 320 Myr in that particular model, while the value to the bottom-left indicates the Ho- emission (mR) at the 90 percentile level. A circle (•) 
indicates whether the maximum Her intensity exceeds 160 mR; a cross (x) signals failure to do so. The circle diameter is roughly proportional to the maximum 
Her intensity in each case. The grey hatched area highlights models that are inconsistent with the halo mass constraint (see text for details). Note that the top 
(bottom) panels correspond to models where n(55 kpc) = 2 X 10“^ cm“^ (4 X 10“^ cm“^), and left (right) panels to models where a low (high) resolution beam 
with diameter 9 = 1° (10’) has been adopted to produce virtual observations. All models at > 100 kpc have been omitted since their respective Her emission 
is weak. 


model. For this particular snapshot and model, we also ob¬ 
tain the value of the emission at the 90 percentile level of the 
corresponding distribution. We then assess the performance 
of each model simply by comparing both the maximum Ha 
emission to the mean level of Ha emission (~ 160 mR) ob¬ 
served along the Magellanic Stream. The result of these ap¬ 
proach applied to both the standard and the extended models 
is collected in Figure 8. 

There, the top panels correspond to the standard models, 
and the bottom panels show the results for the extended set; 
the left (right) panels correspond to the low (high) resolution 
cases. In any panel, each ‘data point’ corresponds to a partic¬ 
ular model identified by {«(55 kpc), 6, t, r„s), and it consists 
of a symbol (circle or cross), and two numerical values. The 
value to the top-right of a given symbol (in parentheses) in¬ 
dicates the maximum Ha emission (in mR), while the value 
to the bottom-left indicates the value at the 90 percentile level. 
A circle indicates whether the maximum Ha intensity exceeds 
160 mR; a cross signals failure to do so. The circle diameter 
is roughly proportional to the maximum Ha intensity in each 


case. In addition, we have greyed hatched the parameter space 
corresponding to models that are deemed ’unphysicaT as per 
the above discussion. As a guide, note that the series of num¬ 
bers in parentheses, i.e. the maximum Ha intensity, shown 
on the top-left panel for model t - I correspond to the results 
shown in the top-right panel of Figure 5 at 220 Myr (ignoring 
the 150 kpc series). Similarly, the maximum Ha intensities at 
r„s= 55 kpc for all models correspond to the results shown in 
the bottom-right panel of Figure 5 at 170 Myr r = 0.5, and at 
220 Myr for all the other models. 

Apparently, both a higher density and an improved resolu¬ 
tion enhance the Ha emission, but in different ways. On the 
one hand, increasing the density shifts the overall Ha intens¬ 
ity towards the high-end. This can be seen by comparing the 
maximum Ha intensity (and the 90 percentile) between the 
standard and the extended models, which are roughly a factor 
2-3 higher in the latter. Increasing the resolution, on the other 
hand, boosts only the maximum Ha emission, without signi¬ 
ficantly affecting the distribution of intensities as a whole, as 
can be judged by comparing the values at the 90 percentile 
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level. 

The standard models at low resolution (top-left panel) fail 
dramatically at any distance in matching the Stream’s mean 
emission. Both the standard model at high resolution (top- 
right) and the extended model at low resolution (bottom-left) 
result in Ha intensities in the 100 mR regime only in the near 
field at 55 kpc. In this sense, both a higher halo density model, 
and virtual observations with a high resolution are equally cru¬ 
cial factors in pushing the Ha emission towards higher levels. 
It is worth mentioning at this point that the original model 
by Bland-Hawthorn et al. (2007) - an instance of a ‘standard’ 
model in our terminology - implicitly assumed an infinite res¬ 
olution, and was therefore capable of reproducing Ha emis¬ 
sion at levels of a few hundred milli-Rayleigh at 55 kpc. 

The extended, high resolution models (bottom-right) are the 
most promising of all the models considered here. In the near 
field, the shock cascade in all these models results in Ha in¬ 
tensities which reach, or even exceed, the highest levels of 
~ 700 mR observed along the Magellanic Stream over the 
SGR In the far field (r„s ~ 75 - 100 kpc), two of the mod¬ 
els (t = 0.5 and 0.75) produce Ha emission consistent with 
the mean emission of ~ 160 mR observed along the Stream. 
However, the success of these models comes at a cost. The in¬ 
crease in the halo density is accompanied by an increase in the 
halo mass. This makes all but the t = 1 model be inconsistent 
with the limit on the gas mass of the Galaxy imposed by the 
cosmic fraction of baryons relative to the total mass. Given 
that the isothermal sphere model yields an upper limit on the 
halo gas density at any distance and for a given gas mass, the 
situation is even more unfortunate for any other reasonable 
DM halo model. 

5. DISCUSSION 

Within the context of the ‘shock cascade’ model, we have 
shown that the interpretation of the MS optical emission (a; 
100 - 200 mR) away from the Magellanic Clouds may still 
work for the updated parameters of the Galaxy under a narrow 
set of conditions. Conventionally, the MS was assumed to 
be on a circular orbit at the midpoint of the LMC and SMC 
{d » 55 kpc). At this distance, for a smooth halo density 
profile, as long as the coronal halo density satisfies 2 x 10“^ < 
(«/ cm“^) < 4 X 10“'^, the shock cascade generates sufficient 
Ha emission to explain the observations. The upper limit on 
density ensures that the mass of the corona (when including 
the other baryonic components) does not exceed the baryonic 
mass budget of the Galaxy. 

In recent years, the first accurate measurements of the MC’s 
proper motions, combined with a smaller estimate of the 
Galaxy’s total mass, have led to a major revision of their bin¬ 
ary orbit about the Galaxy. A highly elliptic orbit is now fa¬ 
voured by most researchers, which pushes the Stream’s mean 
distance further out than the conventional assumption. For a 
smooth halo, the lower density limit of the above range can 
occur at » 75 kpc without violating the constraint on the ba¬ 
ryonic mass budget. If the Stream’s mean distance (especially 
over the SGP) happens to exceed this limit, then we are forced 
to either reject the model, or consider more complex density 
distributions for the Galactic corona (see below). 

An alternative interpretation of the Stream’s optical emis¬ 
sion has recently been put forward in the context of the 
Galaxy’s nuclear activity. The energetic bubbles observed 
with the Fermi-LAT in gamma rays (Su et al. 2010) indicate 
that a powerful event has taken place at the nucleus in the re¬ 
cent past. If the gamma rays are produced through inverse 


Compton upscattering of soft photons, this event can be dated 
to 1-3 Myr ago (Guo & Mathews 2012). Within the context 
of this model, Bland-Hawthorn et al. (2013) show that the 
Stream Ha emission can also be explained by accretion-disk 
driven ionization for Stream distances of 100 kpc or more 
over the poles. The recent discovery of high ionization spe¬ 
cies (e.g. Si IV, Civ) over the SGP (Fox et al. 2015) may 
lend further support to this model. On the other hand, Bland- 
Hawthorn et al. (2007) provide diagnostics of slow shocks 
(e.g. Balmer decrement) that are likely to be observable along 
the Magellanic Stream in future observing campaigns. If en¬ 
hanced Balmer decrements (Ha/HjS > 3) are confirmed along 
the Stream, then some variant of the shock cascade model may 
be needed. 

There is now increasing evidence that the CGM of low red- 
shift galaxies is multi-phase, with a comparable fraction of ba¬ 
ryons both in a hot and a warm phase (Werk et al. 2014). Mod¬ 
ern simulations of the CGM also suggest that the hot halos 
of galaxies are likely to be heavily structured, at least during 
a major phase of gas accretion. While CDM accretion may 
be isotropic on average, individual events involving massive 
systems are not, as clearly demonstrated by the Magellanic 
Clouds in orbit around the Galaxy, or the Sgr dwarf which 
extends through much of the halo (Ibata et al. 1994). The 
mass of this system was probably comparable to the LMC 
and may well have retained gas before being tidally disrupted. 
A more recent accretion event is attested by the massive H i 
stream, the Smith Cloud (> 2 x IO^Mq), that is presently be¬ 
ing stripped and ablated by the corona (Bland-Hawthorn et al. 
1998; Lockman et al. 2008; Nichols et al. 2014). 

We may need therefore to consider the possibility that the 
corona is inhomogeneous rather than smooth. This would al¬ 
low for significant density variations along different directions 
and at different distances, without violating constraints on the 
total baryonic mass of the Galaxy. In turn, based on our exten¬ 
ded models, this would restore the shock cascade as a viable 
model to explain the Ha emission, allowing the Stream to lie 
at ~ 100 kpc the SGP as predicted by dynamical models (e.g. 
Guglielmo et al. 2014). In this scenario, we may envisage 
the strong optical emission along the Stream as a result of the 
Stream’s gas colliding with high-density debris of past accre¬ 
tion events scattered along its orbit. At this point, however, it 
is not clear what a suitable model for an inhomogeneous hot 
halo might be. We will address this in future work. 

If confirmed, the larger distance to the Stream of 100 kpc 
would lift its mass to roughly 8x10® Mq, comparable to the 
total coronal gas mass. If it all breaks down, it would roughly 
double the mass of the corona, at least for a while. We thus 
speculate that the halo of the Galaxy is substantially mass- 
loaded with gas lost by smaller accreted systems. The inter¬ 
action with the hot halo may prevent this gas from cooling 
sufficiently to condense and ‘rain’ down on the disc. Such a 
process is analogous to the meteorologic phenomenon know 
as virga, a type of atmospheric precipitation that evaporates 
while dropping and thus fails to reach the ground. The heavy 
halo may thus serve as a huge reservoir, from which gas may 
eventually be forced out by the strong interaction at the disc- 
halo interface (e.g. Marinacci et al. 2010). 
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APPENDIX 


Ha RECOMBINATION COEEEICIENT 

We describe the temperature dependence of the hydrogen total case B recombination coefficient, as, and of the effective Her 
recombination coefficient, cr^“\ with the generic fitting formula (Pequignot et al. 1991): 


a{T) — ao X 


{l+c)(TI\0^Kf 
\+c{TI\0^KY ■ 


(Al) 


The parameter values appropriate in each case are, respectively, cro = 2.585 X 10“^^ cm^ s“^ b - -0.6166, c — 0.6703, and 
d - 0.5300, and tro = 1.169 x 10“'^ cm^ s“', b - -0.648, c = 1.315, and d - 0.523. Note that this formula is accurate to two 
percent in the range 40 K < T < 2 x 10"^ K. 
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T. Tepper-Garcfa, J. Bland-Hawthorn, and R. S. Sutherland 
PHOTOIONISATION-INDUCED Ha EMISSION 


We include in our model the contribution of the cosmic UV background (UVB) ionising radiation, which leads to a low, but 
non-negligible level of Ha emission along the Stream. Eor simplicity, we assume that the gas is highly ionised («h ~ «Hu); and 
photoionisation equilibrium, which implies that the ionisation and recombination events balance each other: 

(Ue HHii) = HHi Ehi, (B1) 


Here, Thi is the Hi photoionisation rate in units of photon per atom per second; Jy is the angle-averaged specific intensity of 
the UVB; ~ 3.29 x 10'^ Hz is the minimum photon frequency required to ionise hydrogen; and crH(v) a; cro (v/vll)^ is the 
hydrogen ionisation cross section with ctq = 6.3 x 10“'* cm^. 

The Ha emission induced by the metagalactic ionising radiation field along the sightline is thus given by equation (14), where 
the fraction of recombinations that produce an Ha photon is 

fHa(T) = ^ ^ 0.452 g(T ), (B2) 

asiT) 


Here, o' is a monotonically decreasing function of temperature (see equation Al) 


g{T)^ 


/I+c\ 

I H-c' {TIlO'^Kf 

U+c'i 

1 +c{Tll0^K)‘^ 




(B3) 


which satisfies ^(lO^ K) = 1, and g e (0.6, 1.3) for T e [10*, 10®] K. 

In general, /h„ (through T) and «hi both vary along the sightline. However, an estimate of the Ha signal resulting from the 
ionisation by the cosmic UVB of gas at can nevertheless be obtained assuming the gas temperature to be uniform along the 
sightline. In this case, and inserting the appropriate numerical values we get 


~ 452 mR 


Nh. \ 
10** cm“* / 


Th. 

10-'2 S-' 


(B4) 


where Nhi is the integral of «hi along the sightline. Hence, for an Hi ionisation rate Thi ~ 10“'* s“* at z = 0 and gas at 10^ K 
with N Hi ~ 10** cm“ , the Ha signal resulting from the ionisation by the cosmic UVB of gas at 10"* K is roughly 5 mR. 

Since we are not performing proper radiative transfer calculations, we limit the depth (along the sightline) of the gas ionised 
by the UVB to a value Lmax defined by the condition that the column recombination equals the incident ionising photon flux i/r,: 

r*""“ 

as I {ne uhii) ds - ipi , (B5) 

Jo 

where the ionising photon flux (in photons cm“*s“*) is given by equation (15). This condition implies that all the ionising photons 
be absorbed within a depth L^ax, assuming the gas has been exposed to the (uniform) UV radiation field long enough to reach 
ionisation equilibrium (which is well justified in the case of the Stream; see Bland-Hawthorn et al. 2013, their Appendix). Note 
that this condition is equivalent to restricting the ionising effect of the UVB to a column of gas ~ 10** cm“*. Indeed, using 
equations (BI) and (15), the condition (B5) becomes (using y - 1.8) 

2 

iiHids = -Nhi(LL). 


where Nhi(LL) = cTq* » 1.6 x 10** cm *. 









